Retrieve local ijk's of an LGR, given its name #828
Conversation
|
jenkins build this please |
aritorto
force-pushed
the
elemIdxToLocalIJK
branch
2 times, most recently
from
1a44eec to
482e329
Compare
opm/grid/cpgrid/CpGrid.cpp
Outdated
|
|
||
| for (const auto& element : elements(levelGridView(level))) { | ||
| std::array<int, 3> local_ijk; | ||
| levelCartMapp.cartesianCoordinate(element.getEquivLevelElem().index(), local_ijk, element.level()); |
There was a problem hiding this comment.
This looks a bit to complicated. element should be on the level and should be passed directly.
Why do we need to specify the last parameter? Shouldn't the level should be implied by the element.
There was a problem hiding this comment.
Indeed, it's enough with element.index(). The last parameter we need since the LevelCartesianIndexMapper takes a CpGrid in it's constructor, instead of a CpGridaData. Then, you should specify the level.
I'll check any improvements on LevelCartesianIndexMapper to avoid this last argument.
opm/grid/cpgrid/CpGrid.cpp
Outdated
| @@ -713,6 +714,30 @@ std::vector<std::array<int,2>> CpGrid::mapLeafIndexSetToLocalCartesianIndexSets( | |||
| return leafIdx_to_localCartesianIdxSets; | |||
| } | |||
|
|
|||
| std::vector<std::array<int,3>> CpGrid::localIJK(const std::string& lgr_name) const | |||
There was a problem hiding this comment.
Should we make this a free function instead of a member function?
There was a problem hiding this comment.
Yes, we should. Then lgr_names_ could be replace by a call of getLgrNameToLevel()
|
@blattms At some point, some methods related to refinement can be moved into CpGridUtilities, in case the header is kept |
aritorto
force-pushed
the
elemIdxToLocalIJK
branch
from
23d1b90 to
57114b2
Compare
|
jenkins build this please |
1 similar comment
|
jenkins build this please |
…ne - correct localIJK
aritorto
force-pushed
the
elemIdxToLocalIJK
branch
from
2e506d3 to
1db6f0a
Compare
… lgr with inac cell
aritorto
force-pushed
the
elemIdxToLocalIJK
branch
from
1db6f0a to
bd231ba
Compare
|
jenkins build this please |
To facilitate output file generation, it is beneficial to retrieve the local IJK coordinates from each refined level grid based on its name.
Currently, inactive cells in LGRs are not supported in "opm-common." Nevertheless, an lgrIJK() test for a grid containing inactive cells within the LGR block has been added (and commented out). This case test passes if initLgr() is not invoked in opm-common/opm/input/eclipse/EclipseState/EclipseState.cpp.